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Abstract 

We report detailed shape measurements of the tips of three-dimensional NH4CI 
dendrites grown from supersaturated aqueous solution. For growth at small super- 
saturation, we compare two models for the tip shape: parabolic with a fourth-order 
correction and a power law. Neither is ideal, but the fourth-order fit appears to 
provide the most robust description of both the tip shape and position for this ma- 
terial. For that fit, the magnitude of the fourth-order coefficient is about half of the 
theoretically expected value. 
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1 Introduction 



Dendrites are a commonly observed microstructural form resulting from the 
diffusion-limited solidification of non-faceting materials, and they continue 
to be interesting for both practical and aesthetic reasons. Practically, an 
improved understanding of dendritic microstructures may enhance the abil- 
ity to predict and control material properties. Aesthetically, they are an in- 
triguing example of pattern formation under non-linear and non-equilibrium 
conditions. [1,2,3,4] 

Two of the most basic experimental characterizations of a growing dendrite 
are the tip size and growth speed. Although there is already considerable data 
available, recent advances in both theory and experimental technique have 
made more precise comparisons between them possible. Some comparisons 
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Fig. 1. Dendrite of NH4CI. The picture is approximately 400 /jm across. 

have already been made for some pure materials [5, 6], but it is important that 
they also be made for as broad a range of systems as is reasonable, in order 
to clarify the roles of various effects such as solution vs. thermal growth, and 
different values of crystalline anisotropy. 

In the absence of surface tension, one solution to the diffusive growth prob- 
lem is a parabolic dendrite with radius of curvature p propagating at con- 
stant speed v. However, the presence of surface tension, and the instabili- 
ties that lead to sidebranching, complicate the problem considerably. Indeed 
even the most basic issues about the precise tip shape and whether a den- 
drite actually grows with a constant velocity are still areas of active research 
[5,6,7,8,9,10,11,12,13,14]. 

Addressing those issues requires unambiguous ways to identify both the tip size 
and position. In this paper, we consider two different models for the dendrite 
tip shape and evaluate their use to characterize the growth of NH4CI dendrites 
at small supersaturations. A similar exploration of different tip shape models 
for dendrites resulting from phase-field simulations was reported by Karma, 
Lee, and Plapp[ll], but it is important to investigate how well the different 
models work for different materials under actual observation conditions. 



2 Background 

A typical dendrite of NH 4 C1 grown in this study is shown in Fig. 1. Since 
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Fig. 2. Dendrite of NH4CI with <p « 45°. The scale is the same as in Fig 1. 

NH4CI has cubic symmetry, four sets of sidebranches grow approximately per- 
pendicular to the main dendrite stem. In Fig. 1, two sets of sidebranches are 
visible in the plane of the image; two additional branches are growing perpen- 
dicular to the plane of the image. 

The coordinate system used in this work is defined as follows. The main den- 
drite stem is used to define the z axis. The growth direction is taken as the 
negative- z direction, so the solid crystal lies along the positive z axis. The x 
axis is defined as the direction in the plane of the image perpendicular to z. 
Lastly, is defined to be the rotation angle of the crystal around the z-axis. 
A dendrite with 45° is shown in Fig. 2. 

2.1 Parabolic fit with fourth-order correction 

For diffusion-controlled growth in the absence of surface tension, the Ivantsov 
solution is a paraboloid of revolution of radius p growing at speed v. Once 
anisotropic surface tension is included, microscopic solvability [3, 15] predicts 
both the tip size p and speed v depend on the crystalline anisotropy 64. For a 
fourfold-symmetric crystal such as NH4CI, that anisotropy can be expressed 
in spherical coordinates as [12] 

7(0,$)/7o = 1 +e 4 [4cos 4 + (3 + cos 4$) sin 4 0] (1) 

where 7 is the surface free energy, and and $ are the usual spherical angles. 

Although the presence of anisotropic surface tension is critical for the develop- 
ment of the dendritic structure, the overall magnitude of the surface tension 
is small, so the deviations of the tip shape from parabolic might be expected 
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to be rather small as well. Many experiments, including the landmark ex- 
periments of Huang and Glicksman[16], confirm that a parabola is indeed a 
reasonably good approximation to the tip shape. Experiments on NH 4 Br also 
showed that a parabola is a reasonably good approximation, at least relatively 
close to the tip [17]. 

In the limit of small fourfold-anisotropy, Ben Amar and Brener[18] found that 
the lowest order correction to the parabolic shape is a fourth-order term pro- 
portional to cos 40, where is the rotation angle around the z axis. Thus, at 
least close to the tip, the tip shape could be reasonably well described by 



where (x tip , z tip ) is the location of the tip, p is the radius of curvature at the 
tip, and A 4 = 1/96, independent of anisotropy strength[19]. Using a somewhat 
different approach, McFadden, Coriell, and Sekerka[12] found that, to second 
order in e 4 , A 4 = e 4 + 12e 4 , at least very close to the tip. Using an estimate of 
e 4 0.016 (the value reported for NH 4 Br[17]) this corresponds to At ps 0.019. 

Tip shapes consistent with this model were found by LaCombe and coworkers. [5,8] 
They studied succinonitrile dendrite tips under a variety of three-dimensional 
crystal orientations, evaluated both second order and fourth-order polynomial 
fits, and concluded that the fourth-order fit worked significantly better. [5,8] 

2.2 Power-law 

Further back from the tip, the crystalline anisotropy tends to concentrate 
material into four "fins" such that the shape is no longer well-described by 
Eq. (2). The width of the fins is predicted[20] to scale as (z — z tip ) 3 ^. Scaling 
consistent with this prediction was observed in the average shapes of NH 4 C1 
and pivalic acid dendrites grown from solution [2 1,22]. 

Although this power-law scaling was originally proposed to describe the shape 
of the crystal in the region behind the tip, Bisang and Bilgram found that for 
xenon dendrites with ps 0, this power law was a good fit even quite close to 
the tip. [14,23]. Hence the power law offers another way to characterize the tip 
size and location. 

In this model, we describe the tip shape by 




— A 4 cos 40 
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where we have included the factor of 2 by analogy with Eq. 2. The parameter 
p still sets a length scale for the dendritic structure, but the curvature at the 
tip is no longer defined. 

2.3 Experimental Considerations and Model Limitations 

Each of these models has different limitations. From the theoretical perspec- 
tive, the fourth order fit is only appropriate very close to the tip, well before 
sidebranches become significant. Hence in order to avoid contamination from 
sidebranches, only data with (z — z t i P ) < z max should be used in the fit. Since 
sidebranching activity is detectable even close to the tip (at least for NH 4 Br 
dendrites [24]), z max should not be made too large. On the other hand, z max 
should not be made too small since the the region around the tip contains 
the sharpest curvatures and the largest concentration gradients, and hence 
is the the most subject to optical distortions[5,17]. Considering both issues, 
Dougherty and Gollub[17] suggested z max = 3p as a compromise for parabolic 
fits. 

In contrast, the power-law fit is only appropriate further behind the tip, so its 
usefulness for describing the tip size and position must be explicitly tested. 
However, since it is not necessarily constrained to as small z max , the fit can 
include data points less contaminated by optical distortions near the tip. On 
balance, for xenon dendrites with 0^0, Bisang and Bilgram found that this 
power law provided a reasonable fit. 

Both models are potentially sensitive to the choice of z max used in the fitting 
procedure, though such dependence ought to be minimal if an appropriate 
fitting function is used. Singer and Bilgram discuss a procedure to determine 
p from polynomial fits in a way that is somewhat less model dependent [25], 
but that approach did not offer any significant advantage for this system. 

Thus both the fourth-order and the power-law fit provide reasonable fits, at 
least in some cases, but a direct comparison of the two models for the same 
material is required for an accurate assessment. 



3 Experiments 

The experiments were performed with aqueous solutions of ammonium chlo- 
ride with approximately 36% NH4CI by weight. The saturation temperature 
was approximately 65°C The solution was placed in a 45 x 12.5 x 2 mm 
glass cell and sealed with a teflon stopper. The cell was mounted in a massive 
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temperature-controlled copper block surrounded by an insulated temperature- 
controlled aluminum block, and placed on the stage of an Olympus BH-2 
microscope. The entire microscope was enclosed in an insulating box. 

The temperature of the outer aluminum block was controlled by an Omega 
CN-9000 controller to approximately ±1°C. The temperature of the inner 
copper block was controlled directly by computer. A thermistor in the block 
was connected via a Kiethley 2000 digital multi-meter to the computer, where 
the resistance was converted to temperature. The computer then controlled 
the heater power supply using a software version of a proportional-integral 
controller. This allowed very flexible control over not just the temperature, 
but also over any changes in the temperature, such as those used to initiate 
growth. The temperature of the sample was stable to within approximately 
±5 x 10- 4 °C. 

A charged coupled device (CCD) camera was attached to the microscope and 
images were acquired directly into the computer with a Data Translation 
DT3155 frame grabber with a resolution of 640 x 480 pixels. The ultimate 
resolution of the images was 0.63 ± 0.01/xm/pixel. As a backup, images were 
also recorded onto video tape for later use. 

To obtain crystals, the solution was heated to dissolve all the NH 4 C1, stirred 
to eliminate concentration gradients, and then cooled to initiate growth. Typ- 
ically, many crystals would nucleate. An automated process was set up to ac- 
quire images and then slowly adjust the temperature until all but the largest 
crystal had dissolved. 

This isolated crystal was allowed to stabilize for several days. The temperature 
was then reduced by 0.77 °C and the crystal was allowed to grow. The crys- 
tal was initially approximately spherical, but due to the cubic symmetry of 
NH4CI, six dendrite tips would begin to grow. The tip with the most favorable 
orientation was followed, and images were recorded at regular intervals. 

The interface position was determined in the same manner as in Ref. [17]. The 
intensity in the image was measured on a line perpendicular to the interface. 
Over the range of a few pixels, the intensity changes rapidly from bright to 
dark. Deeper inside the crystal, the intensity begins to rise again. (This corre- 
sponds to the brighter areas inside the crystal in Fig. 1.) In the bright-to-dark 
transition region, we fit a straight line to the intensity profile. We define the 
interface as the location where the fitted intensity is the average of the high 
value outside the crystal and the low value just inside the crystal. 

This fitting procedure interpolates intensity values and allows a reproducible 
measure of the interface to better than one pixel resolution. It is also insensi- 
tive to absolute light intensity levels, to variations in intensity across a single 
image, and to variations in intensity inside the crystal well away from the 
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interface. For the simple shapes considered here, this method is more robust 
and requires less manual intervention than a more general contour-extraction 
method, such as the one described by Singer and Bilgram for more complex 
crystal shapes [26]. 

This method works best if the image is scanned along a line perpendicular to 
the interface. Since the position and orientation of the interface are originally 
unknown, an iterative procedure is used until subsequent iterations make no 
significant change to the fit. Specifically, we start with an initial estimate of 
the size, location, and orientation of the tip, and use that to scan the image 
to obtain a list of interface positions up to a distance z max back from the 
tip (where z max is some multiple of p). We next rotate the data by an angle 
9 in the x — z plane and perform a non-linear regression on Eq. 2 or 3 to 
determine the best fit values for x t i P , z tip , p, and At (if applicable), and the 
corresponding x 2 value. We then repeat with different 9 values and minimize 
X 2 using Brent's algorithm. [27] This completes one iteration of the fitting 
procedure. We use this result to rescan the image along lines perpendicular to 
the interface to obtain a better estimate of the interface location and begin 
the next iteration. The procedure usually converges fairly rapidly. Even for a 
relatively poor initial estimate, it typically takes fewer than 20 iterations. 

There are some subtleties to the procedure worth noting. First, for a typical 
crystal in this work such as in Fig. 1, only about 120 data points are involved 
in the fit for z max = 6p. For the fourth-order fit with five free parameters, 
there are often a number of closely-spaced local minima in the x 2 surface, 
with tip positions and radii varying by a few hundredths of a micrometer. If 
the iterative procedure enters a limit cycle instead of settling down to a single 
final value, we select the element from that limit cycle with the minimum 
X 2 ■ Second, it is worth noting that a generic fourth-order polynomial fit is 
inappropriate for Eq. 2, since (after rotation) there are only four parameters 
to fit: Xup, ytip, Pi and A4. Finally, we have no way to control or precisely 
measure the orientation angle of the crystals in our system. 



4 Results 

The best estimates of p as a function of z max /p for the crystal in Fig. 1 are 
shown in Fig. 3. We have included results for the fourth-order and power-law 
fits as well as for a simple parabola for comparison. The corresponding x 2 
values are shown in Fig. 4. 

None of the fits is robust very close to the tip, indicating that the actual tip 
shape is not well-described by any of the candidate functions. The parabolic 
fit also gets rapidly worse for z max greater than about 5 p. The fourth-order fit 
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Fig. 3. Best fit value of p as a function of z max /p for the crystal shown in Fig. 1. 
The curves are for a parabolic fit (+), a fourth-order fit (x), and a power-law fit 

(*)■ 




5 10 15 20 25 

Zrnax/P 



Fig. 4. Values of x 2 for the best fit as a function of z max /p for the crystal shown in 
Fig. 1. The curves are for a parabolic fit (+), a fourth-order fit (x), and a power-law 
fit (*). 
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Fig. 5. Best-fit values for the x-position of the tip as a function of z max / p for the 
crystal shown in Fig. 1. Symbols are as in the previous figures. 

appears to have a plateau between roughly 5 and 8p, but p gradually increases 
with z max , and the x 2 value rapidly increases for z max > lOp. By contrast, the 
power law fit gives relatively stable values at large z max for both p and x 2 ■ 

A second important consideration is the degree to which each fit accurately 
describes the tip location. This is illustrated in Fig. 5, which shows xu p for each 
fit. Here again, the fourth-order fit is slightly more robust than the power-law. 

Finally, Fig. 6 shows the original data along with each of the three model fits 
superposed for z max = 6p. Of the three fits, the fourth-order does the best job 
of capturing both the tip location and shape. 

For the value of z max = 6p, we find A 4 = 0.004 ±0.001, which is similar to that 
measured by LaCombe et al. for succinonitrile[5,8], and to that obtained in 
the simulations by Karma, Lee, and Plappfll]. This value for is less than 
the value of 1/96 predicted by Brener[19], and also significantly less than the 
value of 0.019 estimated by McFadden, Coriell, and Sekerka[12]. It is worth 
noting, however, that these predictions are only intended to be valid close to 
the tip, where our fit is not robust. 

The results are slightly different for the crystal shown in Fig. 2, which has 
(j) pa 45°. The best estimates of p as a function of z max /p for the three fits are 
shown in Fig. 7, and the corresponding x 2 values are shown in Fig. 8. 
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Fig. 6. Best-fit parabolic (solid), fourth-order (dashed), and power-law (dotted) 
curves for z max = 6p, along with the original border (+) for the crystal shown in 
Fig. 1. Near the tip, the parabolic fit is too far to the right, while the power law fit 
is too far to the left. The fourth-order fit matches the tip region fairly well. 
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Fig. 7. Best fit value of p as a function of z max / p for the crystal shown in Fig. 2. 
The curves are for a parabolic fit (+), a fourth-order fit (x), and a power-law fit 
(*). 
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Fig. 8. Values of x 2 for the best fit as a function of z max /p for the crystal shown in 
Fig. 2. The curves are for a parabolic fit (+), a fourth-order fit (x), and a power-law 



Both the parabolic and fourth order models fit reasonably well for z max be- 
tween roughly 5 and 8p. Indeed within that range, for the entire run from 
which Fig. 2 was taken, the value for A 4 is 0.0007 ± 0.0009, consistent with 
zero. In contrast, the power law is a poor fit. 

One significant problem with this measurement is that the image in Fig. 2 is a 
projection of the true three-dimensional shape. This is discussed in Ref. [5,8] 
and in considerably more detail in Ref. [11], but, in general, our findings are 
consistent with those of LaCombe et al. [5,8]. 



5 Conclusions 

We have considered two different models for the tip shape: parabolic with a 
fourth-order correction, and power law. For crystals oriented such that <fi ~ 0, 
both give reasonable fits, though the fourth-order fit is slightly better. For 
rotated crystals, such as those in Fig. 2, however, the fourth-order fit is sig- 
nificantly more robust. 

For the crystal in Fig. 1, the coefficient of the fourth-order term is A 4 = 
0.004 ± 0.001, significantly less than the theoretically expected value. These 
findings are consistent with those of LaCombe et al. [5,8] for succinonitrile 
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dendrites. 



By contrast, the power-law fit was reasonably robust for <f> « 0, in agreement 
with the results of Bisang and Bilgram[14,23] for xenon dendrites, but it did 
not work as well for crystals with <fi m 45°. (One important feature of the 
experiments in Refs. [14,23] was the ability to control the viewing angle and 
hence <fi.) 

One other problem with the power law fit is that it does not accurately describe 
the crystal shape near the tip. Accordingly, it may be more difficult to use a 
power law fit to look for the onset of sidebranching or possible tip oscillations. 

Two remaining issues may be relevant for all of the models. First, the ex- 
tent to which optical distortions near the tip affect the results has not been 
addressed. Specifically, since both the concentration gradients and the interfa- 
cial curvature are largest near the tip, the data closest to the tip are the least 
reliable [5, 17]. Second, the extent to which all of these fits are contaminated by 
early sidebranches needs to be investigated. This may be especially important 
in characterizing the emergence of sidebranches as well as in studies of possible 
tip oscillations. 
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